NMR imaging of the soliton lattice profile in the spin-Peierls compound CuGeC>3 
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In the spin-Peierls compound CuGeC>3, the commensurate-incommensurate transition concerning 
the modulation of atomic position and the local spin-polarization is fully monitored at T — by the 
application of an external magnetic field (H) above a threshold value H c ~ 13 Tesla. The solitonic 
profile of the spin-polarization, as well as its absolute magnitude, has been precisely imaged from 
65 Cu NMR lineshapes obtained for h = (H - H c )/H c varying from 0.0015 to 2. This offers a 
unique possibility to test quantitatively the various numerical and analytical methods developed 
to solve a generic Hamiltonian in 1-D physics, namely strongly interacting fermions in presence of 
electron-phonon coupling at arbitrary band filling. 
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The spin-Peierls (SP) transition observed in quantum 
spin 1/2 antiferromagnetic chains [jj] is the magnetic ana- 
logue of the well-known Peierls transition in quasi-one- 
dimensional (ID) metallic chains. A 3D array of chains of 
spins 1/2, in which the position dependent antiferromag- 
netic exchange interaction J between neighboring spins 
provides a coupling between spins and lattice distortion, 
i.e., phonons, undergoes a second order phase transition 
at T — Tsp from a high temperature uniform, param- 
agnetic and gapless phase, to a low temperature phase 
in which the lattice is dimerized. The corresponding col- 
lective singlet ground state is separated by a gap from 
triplet excitations and the cost in elastic energy is com- 
pensated by the gain in magnetic energy (instead of elec- 
tronic kinetic energy in the Peierls transition) . Applying 
the canonical Jordan- Wigner transformation, the Hamil- 
tonian of such a chain placed in an external magnetic 
field H can be mapped onto that of interacting (spinless) 
fermions coupled to the phonons, H playing the role of 
the fermionic chemical potential. Physical realizations of 
SP systems thus offer the unique possibility to check the 
understanding of this fundamental 1-D Hamiltonian for 
arbitrary band filling, and particularly close to half-filling 
(corresponding to H = 0) where commensurability effect 
are expected to play an important role through the un- 
klapp processes gj. Until recently, physical realizations 
of SP system were limited to organic compounds . The 
recent identification of the first inorganic compound of 
this type ||, namely CuGeOa, and the synthesis of large 
and high-purity single crystals, allowing for example in- 
elastic neutron scattering, has therefore raised an impor- 
tant experimental effort Q. In particular, one of the key 
features of the SP systems is the possibility by applying 
a strong enough magnetic field to cancel the gap between 
the singlet ground state and the lowest triplet excitation. 
This induces a transition to a new magnetic phase in 
which the periodicity of the spin-polarization (i.e., the lo- 



cal magnetization) and of the associated lattice deforma- 
tion is incommensurate (IC) with the underlying crystal- 
lographic lattice. This commensurate-incommensurate 
(C-IC) transition takes place by the introduction of a soli- 
ton lattice. Each soliton corresponds to a rapid change of 
phase between two equivalent commensurate (dimerized) 
domains, and is predicted to carry a total spin | ||. In- 
creasing the magnetic field above the threshold value H c , 
the modulation of the lattice distortion and the local spin 
polarization will continuously evolve from a well-defined 
soliton lattice to planewave regime. 

In this letter, we present a detailed study of 65 Cu NMR 
lineshapes obtained in a CuGeC>3 single crystal in the 
field range 13-26 T which gives a direct and detailed ac- 
cess to the real space profile of the spin polarization in 
the IC phase of CuGe03 for h = (H — H c )/H c vary- 
ing from 0.0015 to 2. The spin polarization profiles have 
been successfully parametrized using Jacobi elliptic func- 
tions, which are the generic solutions of the sine-Gordon 
equation, allowing the determination of both the period 
and the shape of the solitons in the whole range of h. 
These results provide a stringent test for analytical and 
numerical solutions of the Heisenberg SP Hamiltonian, 
enabling to investigate the influence of second neighbors 
interactions, 3-D coupling, or phasons. 

In Fig. |l| is shown the evolution of Cu NMR lineshapes 
when H is varied from H c to ~ 2H C . The data are taken 
on the high frequency satellite of 65 Cu isotope in order 
to avoid overlapping of the lines and ensure maximal res- 
olution. All spectra have been taken at T = 4 K and for 
H parallel to the c-axis of the single crystal, i.e., along 
the chains direction. Reducing T to 1.3 K would increase 
the width of the spectra by only 5 % while keeping the 
same lineshape. For values of H close to H c — 13 T, the 
spectra have been recorded at fixed value of H and sweep- 
ing the NMR frequency (Fig. |], bottom panel), and for 
higher values of H at fixed frequency sweeping the mag- 
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netic field (top panel); the latter method is technically 
more convenient for wide lines, and is used whenever the 
modification of the line-shape during the sweep is negli- 
gible. Splitting of the high field edge singularity in these 
spectra is due to slight misalignment of the single crystal. 
Up to the critical field H c , only a narrow symmetric line 
is seen, corresponding to the homogeneous zero spin po- 
larization in the dimerized phase. Just above, this line is 
strongly broadened and a wide distribution due to soli- 
tons appears. With increasing H, the broadened zero 
spin polarization peak is rapidly suppressed, reflecting 
the disappearance of the nearly dimerized regions and 
the approach towards a planewave-like regime. 

Knowing the hyperfine and quadrupolar tensors expe- 
rienced by the copper nuclei in CuGeC>3 || , the resonance 
frequency of a given copper nucleus at site i on the chain 
can be directly and precisely converted into the time av- 
eraged polarization (S z (i)} of the electronic spin located 
on that site. To the spatial distribution of (S z (i)) in the 
IC phase corresponds directly a distribution of local fields 
A cc (S z (i)), where A cc is the hyperfine coupling constant. 
This gives rise to a particular NMR lineshape which is 
nothing but the histogram, or the density distribution of 
(S z (i)) in absolute units of (S z ), thus allowing a determi- 
nation of the absolute value of the order parameter. For 
a symmetric and periodic 1-D modulation of the spin po- 
larization, as expected in the IC phase of a SP system, 
a simple integration of the density distribution directly 
provides the real space profile over half a period . The 
NMR lineshape presented in Fig. [I] have however been 
analyzed somewhat differently, i.e. fitted to the theoreti- 
cally expected profile of the soliton lattice. In the contin- 
uum approximation, the system is described by a phase 
variable which satisfies a sine-Gordon equation (|^] . The 
periodic solutions are given in terms of the Jacobi elliptic 
functions: sn for a lattice distortion and dn (cn) for the 
average (staggered) component of the spin polarization. 
The total spin polarization can then be written as: 
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where i counts the position of the spins in the chain along 
the c direction, £ is the correlation length defined by the 
corresponding sine-Gordon equation, and R is the ratio of 
amplitudes of the staggered and the average component 
of spin polarization, k is the modulus of elliptic functions, 
and W/2 is the amplitude of the staggered component, 
so that W is just equal to the full span {S z ) max — (S z ) min , 
i.e., to the total width of the NMR spectrum. Here the 
frequency or magnetic field axis of the spectra is con- 
verted to the spin polarization using the known hyper- 
fine coupling constant A cc = 4.6 T 0], When the width 
of the spectra is known and their integral normalized to 
unity, the shape of the spectra is determined by only two 
parameters R and k. For each experimental spectrum, 
these have been determined by the least square fit of the 



theoretical shape of the spin density distribution (convo- 
luted with a convenient intrinsic NMR linewidth to allow 
for an experimental resolution) as shown in Fig. ^. 

In addition, the period L can be deduced from the 
NMR line using the fact that it is directly related to the 
average magnetization m = gfiB{S z ) by the relationship 
(S Z )L = 1. Note that this is radically different from usual 
NMR study of incommensurate phases in dielectrics (^] 
or CDW systems jlO| , in which neither the absolute value 
of the order parameter nor the period can be determined. 
The summation of (S z (x)) over a period leads to the spa- 
tial average spin value (S z ) = %W/AK{k)R, where K(k) 
is the complete elliptic integral. In Fig. [| is shown (S z ) 
as a function of the magnetic field, together with 1/L di- 
rectly measured by X-rays, corrected for the anisotropy 
of the Lande tensor g, since 1/L oc 1/g jyj. The agree- 
ment is very good indeed. However, the range of mag- 
netic field available for X-rays was limited to approxi- 
mately 4% above H c , and only the vicinity of the C-IC 
transition could be explored in pure CuGe03 samples. 
(S z ) is found to be proportional to l/ln[8/(H/Hc — 1)], 
a function which is raising very sharply at H c and then 
rapidly becoming proportional to H . This is indeed pre- 
dicted within the mean field theory description |l3| and 
was already found experimentally in the magnetization 
measurements in other spin-Peierls compounds [jp. 

In addition to the period L, the reconstitution of the 
spin polarization profile as shown in Fig. |^ allows the 
determination of the correlation length £ = L/[4kK(k)] . 
As a function of magnetic field £ is somewhat decreasing 
(Fig. |^). In the field range 1.1-2 H c , the values of £/c 
are found equal to 8-7 in agreement with the theoreti- 
cal predictions. As expected, the amplitude of the spin 
polarization within the soliton, i.e. the width W is only 
slightly dependent on magnetic field (Fig. |J). But the 
value of the ratio R between the amplitudes of the stag- 
gered and the averaged components, which is found to 
decrease from 1.4 at H c to 1 at about 2H C , is definitely 
several times smaller than that predicted either by ana- 
lytical H or numerical solutions Jl4],|l5| of corresponding 
model Hamiltonians for an isotropic exchange. Although 
these results are closer to the predictions corresponding 
to the XY model 113] , there is no reason to doubt that 
the CuGeOs system is in the Heisenberg limit. Note also 
that in organic spin-Peierls compounds, the few NMR 
data available in the incommensurate phase [|l6],[l7j have 
been interpreted as evidence for a staggered component 
much larger than the average one, leading to a symmctic 
lineshape broadening. However, the absence of precise 
knowledge of the hyperfine coupling, and the poor def- 
inition of the lineshapes have precluded more detailed 
analysis so far. 

This strong reduction of the staggered component of 
spin polarization as observed by NMR and compared to 
theoretical predictions strongly suggests, if these latter 
are correct, the existence of some averaging process which 
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is not clearly identified yet. Since a static filtering due to 
a competition between on-site and first-neighbors trans- 
ferred hyperfme field can be excluded Q , it sounds nat- 
ural to consider some dynamical process, remembering 
that NMR only measures time averaged values of the 
spin polarization. It was recently suggested by Uhrig et 
al. that this averaging could be due to the zero point 
motion of phasons, which are acoustic-like excitations al- 
ways present in incommensurate structures. What is the 
amplitude of the pinning gap for these excitations, and 
how they should influence the nuclear relaxations times 
and the lineshapes in the IC phase remains to be inves- 
tigated. 

A salient result of the analysis is the magnetic field 
dependence of the shape of the soliton lattice, which is 
described by the modulus k. It turns out that the data 
for the complementary modulus k' follow a power law 
as a function of reduced field and up to 2H C (our high- 
est experimentally accessible value): k' = (1 — k 2 ) 1 / 2 — 
0.56(H/H C — l) 35 , which remains to be compared to the- 
oretical predictions. Our experimental k(H) dependence 
is in clear disagreement to the old prediction based on 
the Hartree-Fock approximation This is not sur- 

prising, as the customary approximation was to take the 
Hartree-Fock parameter site- independent, thus reducing 
the full Heisenberg Hamiltonian to the XY model, which 
yields to the strong modification of the soliton lattice 
profile. This error has been realized only recently dur- 
ing the analysis of the CuGe03 NMR data [jl9| , and true 
site-dependent Hartree-Fock approximation approaches 
indeed the (more) exact solutions, given either by nu- 
merics or by the analytical continuum description using 
bosonization technique. 

In conclusion, we have followed the evolution of the 
soliton lattice as a function of the magnetic field H in 
the incommensurate magnetic phase of pure CuGe03 by 
recording Cu NMR lineshapes in the field range 13-26 
Tesla, i.e. for h = (H — H c )/H c varying from 0.0015 to 
2. Experimental data have been successfully fitted to the 
analytical solutions of the continuum description, provid- 
ing for each value of the field both the period L = 1/ (S z ) 
of the modulation and the width of the solitons £ [E0| . 
This provides a good overall confirmation of the contin- 
uum description of quantum spin 1/2 chains. The ampli- 
tude of the staggered magnetization is however found far 
smaller than predicted by theoretical models, which could 
be due to the zero point motion of phasons. The mag- 
netic field dependence of £ should also be a stringent test 
for numerical or analytical solutions of the spin-Peierls 
Hamiltonian in presence of frustrating second neighbors 
coupling and 3D interactions. 
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FIG. 1. Magnetic field dependence of 65 Cu NMR lineshape 
in a CuGeC>3 single crystal at 4 K for h = (H — H c ) / H c vary- 
ing in the range 0.0015-2. In the bottom panel, the spectra 
have been recorded by sweeping the frequency. The frequency 
scales of these spectra are mutually shifted in order to com- 
pensate for the variations of magnetic field. In the top panel, 
the spectra correspond to field sweeps. 

FIG. 2. Real space reconstitution of the spin polarization 
profile (left) and the corresponding fit to Jacobi elliptic func- 
tions superposed on experimental spectra (right), as explained 
in the text. Left top panel shows one period of soliton lattice 
taken at field only 0.15 % above H c , with soliton peaks well 
separated by "dimerized" regions. High above H c (bottom 
panels) these regions disappear and solitons strongly overlap 
making the profile more sinusoidal. 

FIG. 3. Field dependence of the period L as deduced di- 
rectly from 65 Cu NMR measurements, compared to the data 
available from X-rays scattering [jH]. For convenience, the 
same data have been plotted in a linear (left and top) and a 
logarithmic (right and bottom) scale. 

FIG. 4. Magnetic field dependence of the fit parameters for 
the NMR lineshapes shown in Fig. [j]. Top panel shows the full 
width W (right scale, triangles) of the spectra and the ratio 
R (left scale, circles) of amplitudes of the staggered (cn) and 
the average (dn) component of the spin polarization. Bottom 
panel shows the correlation length £ (right scale, squares) and 
the modulus k' (left scale, circles) of Jacobi elliptic functions 
reflecting the evolution of the lattice from the extreme soliton 
limit to a more sinusoidal profile. 
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